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ABSTRACT 

We present measurements of the two-point galaxy angular correlation function as a 
function of apparent magnitude, color, and morphology. Our analysis utilizes images 
taken using the UH8K CCD mosaic camera on the CFHT. Six 0°5 x 0°5 fields were 
observed for a total of 2 hours each in I and V, resulting in catalogs containing ~ 25000 
galaxies per field. We present new galaxy number counts to limiting magnitudes of 
I = 24.0 and V = 25.0. We divide each passband sample into intervals of width one 
magnitude. Within each magnitude interval, we parameterize the angular correlation 
function by and find a;(0) to be well described by a power-law of index 6 = 0.8. 

We find the amplitude of the correlation function, Ay^, to decrease monotonically with 
increasingly faint apparent magnitude. We compare with predictions which utilize red- 
shift distributions based on deep spectroscopic observations. We conclude that simple 
redshift-dependent models which characterize evolution by means of the e parameter 
inadequately describe the observations. This is because the predictions do not allow 
for the varying mix of morphologies and absolute luminosities (and hence clustering 
strengths) of galaxies sampled at different apparent magnitudes. We find a strong clus¬ 
tering dependence onV — I color. This is because galaxies of extreme color lie at similar 
redshifts and the angular correlation functions for these samples are minimally diluted 
by chance projections. We find extremely red iV — I = 3.0) galaxies (likely early-type 
galaxies at z ~ 1) to have an A^j about 10 times, and extremely blue {V — I = 0.5) 
galaxies (likely local late types) to have an A^j about 15 — 20 times higher than that 
measured for the full field population. 

We then present the first attempt to investigate the redshift evolution of cluster¬ 
ing, utilizing a population of galaxies of the same morphological type and absolute 
luminosity. We study the dependence of oj{9) on redshift for L* early-type galaxies in 
the redshift range 0.2 < z < 0.9. Although uncertainties are large, we find the evo¬ 
lution in the clustering of these galaxies to be consistent with stable clustering [e = 0 
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for a redshift dependence of the spatial correlation function, ^(r), parametererized as 
^(r, z) = (r/r())~'^(l + We find L* early-type galaxies to cluster slightly more 

strongly (physical correlation length rg = 5.25 ± 0.28 h~^ Mpc assuming e = 0) than 
the local full field population. This is in good agreement with the correlation length 
measured by the 2dFGRS for L* early-type galaxies in the local universe. 

Subject headings: cosmology: large-scale structure of universe — cosmology: observa¬ 
tions — galaxies: photometry — galaxies: evolution 


1. INTRODUCTION 

One of the most interesting and important problems in modern astronomy is that of galaxy 
formation and evolution. Traditionally, galaxy distributions have been quantified using correlation 
functions. Two principle approaches, each with its own advantages and drawbacks, have been 
used to measure the two-point function, which measures the excess probability over random of 
finding another galaxy at a separation r from a given galaxy. One option is to compute the spatial 
correlation function, ^(r), or its Fourier transform equivalent P{k), directly using spectroscopic 
redshifts. In recent years a plethora of redshift surveys culminating in the 2-degree Field Galaxy 
Redshift Survey (hereafter 2dFGRS) and Sloan Digital Sky Survey (hereafter SDSS) have been 
providing accurate pictures of the distribution of galaxies in the local universe. The spectroscopic 
approach, however, is limited by the technical difficulties inherent in measuring spectra for many 
galaxies and the unfeasibility of obtaining both a very deep and very wide sample. 

In order to quantify the clustering of faint galaxies beyond spectroscopic limits the angular 
two-point correlation function, uj{9), has been the approach of choice. The power of uj{9) as a 
diagnostic of the galaxy distribution lies in the simplicity of its application; one is required only to 
count pairs of galaxies at given angular separations and normalize the results with respect to the 
number of pairs expected from a random distribution. The drawback is that angular correlation 
function analyses must rely on accurate observations or models of the redshift distribution of faint 
galaxies to invert uj{9), and hence deduce the three-dimensional correlation length, rg. 

Pioneering studies of ijj{9) using photographic plates were carried out by Groth & Peebles 
(1977); Phillipps et al. (1978); Shanks et al. (1980); Maddox et al. (1990); Bernstein et al. (1994) 
and Infante & Pritchet (1995), among others. In the last decade or so the widespread use of 
charge-coupled device (CCD) cameras have permitted studies reaching far deeper limiting magni¬ 
tudes [passbands in parenthesis] e.g., Efstathiou et al. (1991) [U,B,R,I]; Roche et al. (1993) [B,R]; 
Brainerd et al. (1995) [R]; Hudon & Lilly (1996) [R]; Lidman & Peterson (1996) [/]; Villumsen et al. 
(1997) [R]] Woods & Fahlman (1997) [V,R,I]] Brainerd & Smail (1998) [/]; Postman et al. (1998) 
[/]; Cabanac et al. (2000) [V,I]; Fynbo et al. (2000) [R,I]] McCracken et al. (2000) [B,R,I,K]. 
These investigations were able only to target one or two fields up to about 50 arcmin^ in size. 
Several groups attempted to cover larger areas by mosaicing together many separate pointings 
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(Postman et al.; Roche & Bales (1999) [R]) but these observations reached much shallower limiting 
magnitudes. In recent years, the advent of wide-held mosaic cameras on 4-m class telescopes has 
begun to revolutionize the held, permitting studies of unprecedented depth and areal coverage, 
with the corresponding reduction in variance inherent in covering large areas of the sky (Cabanac 
et ah; McCracken et al. (2001) [I]). 

Spatial two-point correlation functions for local, bright, optically selected samples have been 
determined by many authors over the last few decades. Numerous studies have found ^(r) to 
be well described by a power law, ^(r) = (r/ro)”'*', with slope 7 ~ 1.8 and correlation length 
ro ~ 5 h~^ Mpc for r < 15 h~^ Mpc (Davis & Peebles 1983; Loveday et al. 1995; Norberg et al. 
2001). Tq may be interpreted as the correlation amplitude at 1 h~^ Mpc. Thus the value of vq 
provides a measure of the clustering strength of galaxies in the sample, with a larger value implying 
stronger clustering. 

In the last few years, it has become generally accepted that clustering strengths, as measured 
in the local universe, have a dependence on galaxy morphology (lovino et al. 1993; Loveday et al. 
1995; Hermit et al. 1996; Willmer et al. 1998; Norberg et al. 2002). It has been known since at 
least Davis &: Geller (1976) that early-type galaxies (ellipticals and SO’s) cluster more strongly than 
late types i.e., ro(E) > ro(S). Estimates of the ratio of the correlation strengths [ro(E)/ro(S)]^'® 
for the two types range from ~ 1.2 to ~ 5. Claims for a luminosity dependent component to 
galaxy correlation strengths have been more contentious, although results from the most recent 
surveys (Hermit et al. 1996; Lin et al. 1996; Willmer et al. 1998; Guzzo et al. 2000; Norberg 
et al. 2001 ) do indeed seem to indicate that high-luminosity galaxies are more strongly clustered 
than low-luminosity galaxies. A dependence of clustering strength on both intrinsic luminosity 
and morphological type is to be expected if galaxies are biased tracers (see below) of the mass 
distribution in the universe. Early efforts to quantify galaxy evolution as a function of redshift 
utilizing only one passband proved rather limiting. This is because subsets of any given sample 
selected on apparent magnitude contain galaxies of differing absolute luminosities (and differing 
morphological types) at different redshifts, greatly complicating the analysis. Whenever two (or 
more) passbands have been available, color-selection has often been used. Infante &: Pritchet (1995); 
Neuschaefer & Windhorst (1995); Tandy et al. (1996); Roche et al. (1996); Woods & Fahlman 
(1997); Brown et al. (2000); Cabanac et al. (2000); McCracken et al. (2001) and Zehavi et al. 
(2002) all found red galaxies to cluster more strongly than blue galaxies. This is most likely a 
manifestation of the morphological clustering dependence. A more recent innovation has been the 
use of multi-passband data to estimate photometric redshifts (Connolly et al. 1998; Arnouts et al. 
1999a; Brunner et al. 2000; Teplitz et al. 2001), a technique which bridges traditional spectroscopic 
and photometry camps to study galaxies in the range 0 < 2 ; < 1. However, although promising, to 
date, photometric redshifts analyses have been limited to small fields-of-view. 

The “e” formalism, first introduced by Groth & Peebles (1977), has traditionally been used 
to characterize the evolution of clustering with redshift. This empirical approach assumes that 
the typical clustering length observed at high redshift transitions monotonically to that observed 
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locally. By assuming a redshift distribution and cosmology, it is possible to predict compare to 
observations and hence determine the value of e which best describes the evolution of the correlation 
function. Many authors have concluded, either from the aforementioned studies or directly 
from spectroscopic studies (Le Fevre et al. 1996; Shepherd et ah 1997; Small et al. 1999; Carlberg 
et al. 2000; Hogg et al. 2000), that 0 < e < 2 i.e., that galaxy clustering is either stable or 
grows in amplitude from z = 1 to the present. Undoubtedly, a lack of consistency of sample i.e., 
an inability to follow the same population of galaxies due to changing morphological mixes and 
intrinsic luminosities with redshift, makes interpretation of the measurements confusing. 

Strong evidence that the “e” formalism might be invalid, at least at higher redshift, was 
provided by the discovery that clustering amplitudes measured for Lyman-break galaxies at z ~ 3 
are similar to those measured for local galaxies (Adelberger et al. 1998; Giavalisco et al. 1998). Such 
strong clustering of the Lyman-break galaxies is to be expected if these galaxies are highly biased 
with respect to the mass distribution. In the standard hierarchical picture of galaxy formation and 
evolution, galaxy clustering traces overdense regions in the dark matter distribution. High redshift 
galaxies are expected to form at the most extreme peaks in the density field and are thus biased 
tracers of the mass (Kaiser 1984; Bardeen et al. 1986). If early-type or more luminous galaxies are 
associated with rarer, more massive halos, then these galaxies would be expected to exhibit even 
stronger clustering than the galaxy population as a whole. Subsequent to the epoch of formation, 
the clustering of the galaxies is expected to evolve more slowly than the clustering of the dark 
matter so the two distributions are expected to be more similar today than in the past (Baugh 
et al. 1999; Kauffmann et al. 1999; Carlberg et al. 2000). 

In the hierarchical formation scenario, the clustering of dark matter increases monotonically 
with time. The rate of evolution in the clustering of the dark matter is a function of cosmology, 
being faster in a high density (e ~ 1.0 for = 1) than in a low density (e ~ 0.2 for = 
0.2) universe (Cohn et al. 1997). If the evolution in dark matter clustering could be measured 
directly as a function of redshift then flm, ^\, and the power spectrum of density fluctuations 
could be inferred directly. However, in practice, one can measure only evolution in the galaxy 
clustering pattern. The amplitude of the galaxy correlation function is determined by a combination 
of factors: evolution in the underlying dark matter fluctuations plus any bias (which is also a 
function of cosmology) relating the galaxy overdensities to the mass. Realistically, in addition to 
the depth of the dark matter potential, one can expect galaxy evolution also to depend on complex 
physical processes involving the local environment, cooling and feedback mechanisms, and galaxy 
interactions. Observations of galaxy clustering at high redshift are therefore vital to constrain 
models of these processes empirically, and deepen our understanding of the complexities of both 
galaxy and structure formation and evolution. 

In this paper we investigate galaxy clustering evolution on scales of up to 30' using data 
collected using the University of Hawaii’s 8K (UH8K) CCD mosaic camera on the Canada-France- 
Hawaii Telescope (CFHT). The data were obtained for a weak leasing study of “blank fields” i.e., 
the helds chosen for study were intended to be representative views of the universe not containing 
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any unusually large masses such as rich clusters. The data are also, therefore, well suited to the 
study of (field) galaxy clustering and evolution. 

The outline of the paper is as follows. In §2 we describe the data and photometry. In §3 we 
present both our number counts and a comparison with the counts from other groups. In §4 we 
describe and discuss our measurements of the two-point correlation function as a function of median 
apparent magnitude, V — I galaxy color, and morphological type. We compare to predictions and 
discuss possible sources of uncertainty. In §5 we briefly summarize our conclusions. We assume a 
flat lambda (IlmO = 0.3, Hao = 0.7) cosmology with Hq = 100 h km s“^ Mpc“^ throughout. 

2. THE DATA AND GALAXY SAMPLES 
2.1. Data Acquisition and Reduction 

The data were taken at the 3.6m CFHT telescope using the 8192 x 8192 pixel UH8K camera 
at prime focus. The field of view of this camera is ~ 30^ with pixelsize 0.207^h The data (six 
pointings) used in the analysis were acquired as part of an ongoing project that has the principle 
aim of investigating the cosmic shear pattern caused by gravitational leasing from the large-scale 
structure of the universe. In addition to the main project, estimates of cosmic shear variance on 
2' —30' scales (Kaiser, Wilson, & Luppino 2000, Paper I), the data have also been used to investigate 
galaxy halos at radii of 20" —60" (50 — 200 h~^ kpc) (Wilson et al. 2001, Paper II), and to investigate 
the relationship between mass and light on group and cluster scales (45" — 30') (Wilson, Kaiser, & 
Luppino 2001, Paper III). In this paper we focus on properties of galaxy clustering. 

Table 1 gives an overview of the data, describing the field name, center and seeing for each 
pointing. The reduction procedure was lengthy and involved. We defer a full description of the 
data reduction pipeline (involving careful registration and point-spread function correction) and 
the resulting catalogs to a later paper (Wilson & Kaiser 2003). Here we provide only an overview 
of the photometry pipeline. 


2.2. Photometry 

The data were dark subtracted, flat-fielded, registered, median averaged and corrected for 
galactic extinction using extinction measurements from Schlegel, Finkbeiner, & Davis (1998). The 
UH8K science images were calibrated to the Johnson-Cousins system using a series of standard 
star field observations (Landolt 1992). The standard star fields were also dark subtracted and 
flat-fielded in a similar manner to the science frames. These were used then to calculate zero-points 
for each pointing. In agreement with McCracken et al. (2001) our observations did not indicate the 
presence of a color term for either the I or V filters. 
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The imcat data reduction package (http;//www.ifa.hawaii.edu/~kaiser/) was used through¬ 
out. Objects were detected and assigned an optimal size by smoothing the science frames by a 
progressively larger series of filters as described in Kaiser et al. (1995). Aperture magnitudes were 
then measured within a multiple of three times this radius. 

The number counts of objects at faint magnitudes are dominated by galaxies but at brighter 
magnitudes the stellar component is dominant. We removed stars at brighter magnitudes (/ < 23 
; V < 24) by hand, by means of filtering on a size-magnitude diagram. For fainter objects, no 
attempt was made to further eliminate stars from the sample since compact galaxies could be 
mistakenly removed and stellar numbers are very small relative to the galaxies at these faint limits. 
After stellar filtering, approximately 25000 galaxies remained in each passband for each of the six 
pointings, an extremely deep and wide-area dataset compared to previous studies. 


3. NUMBER COUNTS 

The /-band number counts (logarithm of number of galaxies degree”^ mag“^) for each pointing 
are shown in Figure 1. The uncertainties shown are Poissonian. Clearly, the counts from each 
pointing are in good agreement with each other and are complete to / = 24.0. Figure 2 shows the 
same but for the U-band counts which are complete to U = 25.0. 

Figure 3 compares our / number counts with those measured by other groups. Our data are 
shown by the filled hexagons. The measured values may be found in Table 2. The uncertainties at 
each magnitude were calculated from the pointing-to-pointing variations. Also shown in the figure 
are the counts of Postman et al. (1998); Arnouts et al. (1999b); Cabanac et al. (2000); Metcalfe 
et al. (2001); and McCracken et al. (2001). A correction of / = Tab — 0.43 was applied to the 
counts of McCracken et al. to bring them to the Cousins /-band system. Clearly, the counts are in 
reasonably good agreement over the whole range. The highest and lowest values at each magnitude 
were obtained by Cabanac et al. and Postman et al. respectively. However, as discussed by 
Cabanac et al., that group systematically measure 20 — 30% higher counts than Postman et al. but 
have large (0.2-0.5) errors in their zero-point calibration. 

The line in Figure 3 shows the best fit (Table 3) to our counts in the range 21.0 < / < 23.5 
(slope of 0.31 ± 0.01). A comparison between our best-fit slope and the values found by other 
authors is shown in Table 4. Although the range over which the slope is fit varies slightly, the 
values are generally in good agreement. 

Figure 4 shows a similar comparison, but for the (Johnson) V counts (Note that Vab = Johnson 
V). The filled hexagons in Figure 4 represent our measurements. The values are given in Table 5. 
U-band galaxy number counts are more rarely published in the literature. We compare to the 
counts of Gardner et al. (1996); Arnouts et al. (1999b); and Cabanac et al. (2000). We find that 
we are in good agreement with Gardner et al. at bright magnitudes and Arnouts et al. at faint 
magnitudes. However, the agreement with Cabanac et al. is rather poor. Cabanac et al., however, 
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Fig. 3.— Comparison of our /-band number counts with measurements from other groups. The 
filled hexagons represent our data. The uncertainties are obtained from the pointing-to-pointing 
variations in Figure 1. See Table 2 for values. The line shows the best fit to our counts in the range 
21.0 < / < 23.5 (slope of 0.31 ± 0.01). 
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also report large uncertainties in their y-band calibration (ranging from 0.2 at bright magnitudes 
to 0.5 at V of 24). 

The line in Figure 4 shows the best fit to our data in the range 22.0 < V < 24.5 (slope of 
0.43 ± 0.02). Within uncertainties, this value agrees well with the slopes reported by Woods & 
Fahlman (1997); Driver et al. (1994) and Smail et al. (1995) (See Table 4). Our data do not reach 
a sufficiently faint limit to detect any flattening of the slope reported at D ~ 24.5 by Smail et al. 


4. THE ANGULAR CORRELATION FUNCTION 

The two-point angular correlation function, ^(0), measures the excess probability (over a 
random Poisson distribution) that two galaxies will be found in the solid angle elements dfii and 
dfi 2 , separated by angle 6. It is defined by 

dP = N^[i + uj{e)]dnidn2 (i) 

where dP is the joint probability of finding galaxies in the two solid angle elements, and N is the 
mean surface density of galaxies. 


4.1. The Dependence of uj{9) on Apparent Magnitnde and Passband 

We compute uj{9) using the Tandy & Szalay (1993) estimator. 


DD - 2DR + RR 
RR 

where DD, DR and RR are the number of data-data, data-random, and random-random pairs 
(scaled appropriately by the number of data and random points) at angular separations 9 ± d9 
respectively. This estimator is based on the N-point function (Szapudi & Szalay 1998) and has 
been shown to have the advantage of reduced edge effects and smallest possible variance. 



To determine DR and RR we generated a catalog for each pointing containing 50000 random 
points covering a similar area to the data. We masked this random catalog with the same masks 
we used to mask out saturated stars and sub-quality regions of the CCDs. The remaining catalogs 
contained ~ 30000 randomly distributed points. 

We also applied an integral constraint (IC) correction (Groth & Peebles 1977). Estimating the 
mean galaxy density and the two-point correlation function from any survey limited in area results 
in a CO(9) artificially reduced by the amount 


C = 


1 


co{9) dDi dQ .2 


(3) 
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where the integrals are performed over the total solid angle, fl, subtended by each of the pointings 
after masking by the detection masks (Roche &: Bales 1999). The correction is a function of 
the effective survey area and, to a lesser extent, the survey geometry. Because the mean density 
measured from a large CCD more closely approximates that of the mean global value, the correction 
is smaller for larger fields of view. 

In practice uj{6) is calculated without the correction, equation (3) is integrated over all elements 
of solid angle dDj in the survey area and then uj{9) is recalculated with the correction added. A 
stable solution is reached by iteration. In order to calculate the correction it is therefore necessary 
to assume a functional form for the two-point correlation function. For local, bright, optically 
selected galaxy samples lv{9) has been shown to be well approximated by a power law 

aj{9) = A^9-^ (4) 

with 5 = 0.8, and where is the amplitude of uj{9) (Measuring angular separation in arcmin as 
we do here, A^^ is then defined as the amplitude of (jj{ 9) at 9 = 1'). 

For each of our pointings, the IC correction, C, was found to be ~ 0.162A^, varying slightly 
depending on the field geometry and detection mask. The values of C determined for each of the 
six pointings were comparable since the field sizes and geometries were very similar. Note that the 
correction only becomes important when measuring the correlation function for galaxies at large 
separation; for the fields-of-view analyzed here, the correction has a very small effect on A^^. 

To calculate an error estimate we utilize the fact that we have six separate pointings and 
compute the uncertainties from the field-to-held variations at each separation. This method of 
determining the uncertainties is superior to most other analyses which are limited to one or two 
pointings. These are forced to employ bootstrap resampling techniques in order to estimate uncer¬ 
tainties. However, in this analysis, correlations between measurements on different scales are not 
taken into consideration, which may result in an underestimate of the x^- 

Figure 5 shows the two-point angular correlation function for four /-band slices, each one- 
magnitude wide, in the range 20.0 < / < 24.0. We calculated uj{9) using logarithmically spaced 
bins of width Alog0 = 0.2. On small scales (log]^Q(0) < 0.2, 9 < 1.58') we estimated the two-point 
function (equation 2) directly from pair counts. On larger scales, uj{9) was determined from counts 
in cells. The four lines in Figure 5 show the best ht to the data in the range — 1 < log]^Q(0) < 1 
(6" < 0 < 10' ), assuming 5 = 0.8. (We also htted A^^ and the slope 5 separately, both with and 
without an IC correction, and found 5 = 0.8 to be a good fit to the correlation function for each 
of the four magnitude intervals). The best-fit value of log^g the logarithm of the amplitude of 
u){9) at 1', as a function of median I magnitude is shown in Table 6. We find a monotonic decline 
in A^ with increasingly faint median I magnitude. 

Figure 6 shows our measurements of uj{9) but for the H-band, for four one-magnitude wide 
slices in the range 21.0 < V < 25.0. As with the /-band, we find a monotonic decline in the 
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amplitude of the two-point correlation function with increasingly faint V magnitude. The best-fit 
values of logj^Q^o; as a function of median V magnitude may also be found in Table 6. 

Figure 7 compares our measured values of clustering amplitude at 1 arcmin (as a function 
of median Cousins I magnitude) with the values obtained by other groups. We compare with the 
results of Efstathiou et al. (1991); Lidman & Peterson (1996); Woods & Fahlman (1997); Brainerd & 
Small (1998); Postman et al. (1998); Cabanac et al. (2000); McCracken et al. (2000) and McCracken 
et al. (2001). At bright magnitudes (/ < 20) there is good agreement among the measured values 
of Ai^. However, at fainter magnitudes some discrepancies appear. Our measurements are in very 
good agreement with those of Efstathiou et al. (using the value quoted in Lidman & Peterson Table 
6), McCracken et al. (2000) and McCracken et al. (2001). The clustering amplitude we measure 
is stronger than that found by Lidman & Peterson. We measure a somewhat lower clustering 
amplitude than Woods & Eahlman and Cabanac et ah. We measure a significantly lower clustering 
amplitude than Brainerd & Small and Postman et al. 

The discrepancy between our results and those of Woods & Eahlman and Brainerd & Small 
may be explained by the relatively small area investigated by those studies, and the corresponding 
inherently large uncertainties. Woods & Fahlman had field sizes of 3 x 49 arcmin^ and IC correction 
~ 0.02. The corresponding values for Brainerd & Small were 2 x 30 arcmin^ and IC ~ 0.02. 
Cabanac et al. report large uncertainties in their photometric zero points which increase at fainter 
magnitudes. Any variations in effective depth due to zero-point errors or variations in observing 
condition (e.g., variations in the seeing) between exposures would artificially mimic large-scale 
power. The discrepancy with the results of the survey of Postman et al. is more puzzling. The 
large survey of Postman et al. consisted of 16 deg^ (256 x 16 arcmin^, IC ~ 0.002). They find 
a shallower best-fit slope to the correlation function (5 ~ 0.7) for I > 22. This disagreement 
notwithstanding, it seems likely that zero point variations from frame-to-frame may remain in the 
Postman survey. 

Figure 8 compares our measured values of clustering amplitude at 1 arcmin in the E-band with 
the values obtained by Woods & Fahlman (1997) (the only other measurements of published 
in this passband). The uncertainties associated with the Woods Sz Fahlman measurements are 
larger than for our measurements but the agreement is generally good at E < 24. We measure 
a significantly (factor of two) lower clustering amplitude than the Woods & Fahlman faintest 
determination (at E = 24.3). 

In the next section we predict the correlation function amplitude that we would expect to 
measure as a function of apparent magnitude, based on deep spectroscopic observations. We then 
compare these predictions to our measurements. We investigate whether our observations can be 
matched to the predictions assuming simple analytic models for galaxy clustering evolution. 



Fig. 5.— The logarithm of the angular correlation function, to{6), as a function of the lo^ 
of the angular separation in arcmin for various /-band slices as shown. The uncertaini 
calculated from the pointing-to-pointing variations. The lines show the best fits (Table 6( 
range —1 < log^o(^) < 1 < ^ < lOO assuming a power-law slope of —0.8 for uj{6). 
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Median 

Fig. 7.— The logarithm of the amplitude of the angular correlation function uj{9) at 1' for the 
/-passband plotted as a function of median magnitude (see Table 6). The hlled hexagons represent 
our data. Also shown are the measurements from other groups. 
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4.2. Clustering Strength Predictions 

As mentioned in § 1 the two-point spatial correlation function, ^{r), for local, optically selected 
galaxies has been shown to be well described by a power law 


^(r) = (r/ro) (5) 

where r is physical (proper) separation and vq is the correlation length at z = 0. Robust deter¬ 
minations for the power-law slope of 7 ~ 1.8 and physical correlation length of tq ~ 5.0 h~^ Mpc 
have been made by many groups (Davis & Peebles 1983; Loveday et al. 1995; Postman et al. 1998; 
Willmer et al. 1998; Carlberg et al. 2000; Norberg et al. 2001). 

Groth & Peebles (1977) proposed the “e-model” to describe the evolution with redshift of the 
correlation function, measured in terms of proper separation: 


(l + z)-(3+^) (6) 

There are several noteworthy values of the evolutionary parameter e. The clustering pattern is 
fixed in comoving coordinates if e = —1.2 (e = 7 — 3): galaxy clusters expand with the universe 
and clustering does not grow with time. If e = 0.0, “stable clustering”, then clustering is fixed in 
proper coordinates. In this case, the galaxies are dynamically bound and stable at small scales. The 
clustering grows in this case because the background density of galaxies is diluted by the expansion: 
it is effectively the voids that are growing. Linear theory predicts e = 0.8 to describe the evolution 
in the clustering pattern of dark matter in an Dm = 1 universe (Baugh et al. 1999) (See also e.g., 
Colin et al. (1997); Carlberg et al. (2000) and Kauffmann et al. (1999) for e predictions for the dark 
matter evolution in alternative cosmologies). If e > 0 then clustering grows with time in proper 
coordinates, as expected from gravitational collapse. 

The relation then (for small angles) between the two-point angular and spatial coordinates, 
and ^{r,z) is given by Limber’s equation (Limber 1953; Peebles 1980) 


u;{9) 


r(l/2)r[(7 - l)/2] A 
r[(7/2)] 07-1 0 


(7) 


where the Gamma function factor equals 3.68 for 7 = 1.8. A, the amplitude factor, is given by. 



' dN{z)' 
dz 

2 

dz 


roo 

Jo 

dN(z) , ■ 
dz 

■ 2 


( 8 ) 


Here, dN/dz is the number of galaxies per unit redshift interval. The function g(z) depends 



only on 7, e, and the cosmology 
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g{z) = F{x){l + z) 

F gives the correction for curvature, 

F{xf = l + ^k{^^ 

X is the comoving distance at redshift z, 

_ c r 1 

'W) 

E is given by 

1 

E = |^nni(l + z)^ + + z)'^ + ^ 

and rifc (the “curvature of space”) is defined by + + = 1. (Carroll et al. 1992; Hogg 1999; 

Brown et al. 2000) 

Therefore, predictions of galaxy clustering strength, u}{9) are dependent on a number of factors. 
Namely, these are 7, ro, e, the assumed redshift distribution N{z), the Hubble constant Hq, and 
one’s choice of cosmology. If the spatial two-point correlation function is well described by a power 
law (equation 5), then the angular two-point correlation function will also be well described by a 
power law of slope (5 = 7 — 1 (equation 7). The physical correlation length at 2; = 0, ro (to the power 
7), determines the normalization of uj{9). In making our predictions we adopt a local correlation 
length of ro = 4.9 h~^ Mpc, the best-fit real space tq determined by the 2dFGRS (Norberg et al. 
2001 ). 

Accurate predictions of the correlation amplitude at each magnitude depend crucially on the 
assumed redshift distribution (Note that uj{9) decreases with increasing apparent magnitude be¬ 
cause of the increasing probability of chance projected alignments of galaxies at very different 
redshifts). uj{9) is extremely sensitive to the shape of the assumed redshift distribution but not 
to its normalization. As in Wilson et al. (2001), in making our predictions, we adopt redshift 
distributions appropriate to each magnitude interval based on spectra from the SSA22 field sample 
of Cowie (Cowie et al. 1994; Cowie et al. 1996; Cowie et al. 1999; Wilson et al. 2002). We model 
the normalized redshift distribution as 

p{z) = QFz^ ex.Y>{-z/zQ)/zl (13) 

for which the mean redshift is T = Szq and the median redshift is ^median = 2.67zq. Note that 
equation (13) has only one free parameter, the redshift scale parameter zo- 

We calculate zq (in half-magnitude intervals) by setting the mean redshift to match that from 
the Cowie sample (allowing for varying numbers of galaxies in each magnitude interval). The Cowie 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 
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sample reaches a limiting magnitude of / ~ 23.5, V ~ 24.5. At the faintest magnitudes the sample 
becomes about 20% incomplete. It is thought that the galaxies for which a redshift cannot be 
determined lie predominantly around z = 1.5 — 2.0. We therefore calculate zq in two ways. Firstly 
we utilize only galaxies in the Cowie sample with secure redshifts, our “raw” model. Secondly 
we assign a redshift of 1.8 to the galaxies from the Cowie catalog without secure redshifts, our 
“corrected” model. This increases the value of the redshift scale parameter zq we employ, only at 
the faintest magnitudes where incompleteness becomes important. It is likely that the true redshift 
distribution lies between these two extremes. The value of zq we adopt to describe the redshift 
distribution in each half-magnitude interval is shown in Table 7. 

In making these predictions of we assume a flat lambda (fimo = 0.3, Hao = 0.7) cosmology 
with Hq = 100 h km s~^ Mpc“^. The predictions are rather insensitive to the choice of cosmology, 
because the majority of galaxies lie at relatively low redshift, especially at bright magnitudes. 
Adopting an Einstein de-Sitter cosmology would decrease the predictions by about 0.2 in the log 
at the faint end. 

The upper panel of Figure 9 shows the /-band clustering amplitudes from Figure 7 (for clarity 
we do not plot the values obtained by other groups). Also shown are the predictions for assuming 
redshift distributions based on the Cowie sample and also assuming an evolutionary parameter of 
either e = —1.2 (solid), 0.0 (dashed), 1.0 (dot-dashed) or 2.0 (dotted). The upper line in each case 
shows the predictions assuming the raw redshift distribution, the lower line the incompleteness- 
corrected distribution. Clearly, equation (6) does not provide a very satisfactory fit to the data for 
any value of e, assuming a smooth extrapolation of the clustering evolutionary model out to this 
study’s fainter magnitude limits. We find that a clustering pattern fixed in physical coordinates, 
e = 0, provides a reasonable fit to the data at the bright end, but that the data fall below the 
predictions at the faint end. 

The lower panel of Figure 9 shows the predictions for A^ for our E-band data (see Table 6 
for measured values of A^ and Table 7 for the zq adopted to describe the redshift distribution of 
each magnitude interval). As for the /-band, we conclude that equation (6) does not provide a very 
satisfactory fit to the data for any value of e (assuming a smooth extrapolation of the clustering 
evolutionary model out to this study’s fainter magnitude limits). 

In calculating the predicted correlation amplitude we utilized a correlation length of tq = 
4.9 h~^ Mpc as measured by the 2dFGRS. Adopting a larger or smaller local correlation length 
would shift the predictions vertically upwards or downwards. 

We conclude that equation (6) does not provide a good fit to the observed evolution in clustering 
amplitude with redshift. This conclusion was also reached by McCracken et al. (2001). We interpret 
the failure of the predictions to match the data as being due to the assumptions of the model being 
overly simplistic. Equation (6) implicitly assumes that galaxies with different morphologies have 
similar intrinsic clustering properties. This is known not to be the case locally. Davis & Geller 
(1976); Giovanelli et al. (1986); Maurogordato & Lachieze-Rey (1991); lovino et al. (1993); Loveday 
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Median V 

Fig. 9.— Upper panel shows data (solid hexagons) as in Figure 7 for the /-band. The lines show our 
predictions assuming redshift distributions based on deep spectroscopic observations, correlation 
length ro = 4.9 h~^ Mpc as measured by the 2dFGRS, and evolutionary parameter e = —1.2 
(solid), 0.0 (dashed), 1.0 (dot-dashed) and 2.0 (dotted). Clearly, equation (6) does not provide a 
satisfactory fit to the data for any value of e. Lower panel shows same but for the U-band data as 
in Figure 8 (see Table 6 for values). 
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et al. (1995); Hermit et al. (1996); Guzzo et al. (1997); Willmer et al. (1998) and Norberg et al. 
(2002) have found early-type galaxies to cluster more strongly than late types. The morphological 
mix of galaxies in any sample will vary with apparent magnitude due to different K-corrections for 
early and late type systems. As one probes fainter in apparent magnitude, a higher preponderance 
of late types enter the sample (Magliocchetti et al. 2000). Moreover, by selecting on increasingly 
faint apparent magnitude one samples intrinsically fainter and likely more weakly clustered galaxies. 
(McCracken et al. 2001). 

We further investigate the bivariate dependence of galaxy clustering evolution on luminosity 
and morphology in the remainder of this paper. In the next section we investigate uj{9) as a function 
oi V — I color. 


4.3. The Dependence of uj{9) on V — I Color 

In this section we investigate the dependence of clustering strength on galaxy V — I color. 
We employ the same catalogs as used in Wilson et al. (2001), containing galaxies which have been 
detected in both I and V images above a threshold significance of 4zz (to ensure that any given 
“detection” is truly a real object). Galaxies tend to be detected at higher significance in the /-band 
images. We subdivided the data into seven intervals with V — I color ranging from 0.0 to 3.5 and 
A(H — /) = 0.5 (Table 8 shows the number of galaxies in each color interval). As in section 4.1 we 
again assume that oj{9) is well described by a power law of slope 5 = 0.8 and we fit over the range 
— 1.0 < log]^Qa;(0) < 1.0. 

Figure 10 shows clustering amplitude as a function of color for galaxies in the range 20.0 < 
I < 23.0. Figure 2 shows a turnover in the H-band counts at H ~ 25 and thus we may be subject 
to missing fainter galaxies from the redder (V — I >2) intervals. Table 8 shows the best-fit values 
of obtained for each color interval. As in § 4.1, the uncertainties were estimated from the field- 
to-field variations between our six pointings. The dashed line in Figure 10 shows the correlation 
amplitude (—1.80 for a median I magnitude of 22.59) obtained for the full field sample in § 4.1 
(Table 6). 

As seen from Figure 10, red galaxies (V — I c:^2>) have a clustering amplitude about 10 times 
larger than that measured for the full sample in § 4.1. As we shall discuss further in § 4.4, galaxies 
with color V — I = 3 are exclusively early types occupying a narrow range in redshift at about 
z = 1. Since the redshift distribution N(z) for these red galaxies is rather narrow, chance alignments 
of galaxies of the same color at significantly higher or lower redshift cannot occur and hence the 
measured angular correlation function is not diluted as it is for galaxies with colors in the range 
1 < {V — I) < 2, which are more widely distributed in redshift. The high clustering amplitude 
found for extremely red galaxies is also due in small part to these being of a morphological type 
which (at least locally) cluster more strongly than later types. 

As noted by McGracken et al. (2001), we also find an upturn in the clustering amplitude of 
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Fig. 10.— The logarithm of the amplitude of the angular correlation function uj{9) at 1' as a function 
oiV — I color for galaxies in the magnitude range 20 < / < 23 (see Table 8). Extremely red and 
extremely blue galaxies are seen to cluster more strongly than those of intermediate color. The 
dashed line shows the amplitude measured in § 4.1 for the full field sample with median I magnitude 
of 22.59. In all cases, the clustering amplitude is higher than for the full field population. 





- 24 - 


extremely blue galaxies (F — / ~ 0.5). A similar effect was noted by Landy et al. (1996). The 
most likely explanation for the strong angular clustering amplitude found for the blue galaxies is 
that these galaxies are all situated at relatively low redshifts. As with the red sample, the redshift 
distribution N(z) of these extremely blue galaxies is rather narrow and thus is not diluted as 
strongly in projection over redshift as is the intermediate sample. 

The results of this section and of § 4.2 strongly suggest that studying and interpreting galaxy 
clustering either as a function of median magnitude or of color is greatly complicated by the varying 
morphologies and intrinsic luminosities of galaxies contained in the different samples. A far more 
promising approach would be to isolate the same population and study its clustering evolution 
with redshift. Given fluxes in only two passbands, as we have here, it is impossible to select and 
measure evolution in uj{9) for all morphological types. However, as we show in the next section, it 
is possible to use crude photometric redshift determinations to isolate and analyze the clustering 
evolution of L* early-type galaxies with redshift. 


4.4. The Dependence of u}{9) on Morphological Type 

As in Wilson et al. (2001), we use V — I color combined with an /-band magnitude cut to 
select a sample of bright early-type galaxies. This technique depends on the fact that early-type 
galaxies are the reddest galaxies at any given redshift. By selecting galaxies of some color we see a 
superposition of early types at redshift ze such that c = ce{ze) and later types at their appropriate, 
but considerably higher, redshift. An L ^ L-i, early-type galaxy appears much brighter than an 
L ^ spiral galaxy, and as explained in more detail in §2.2 of Wilson et ah, with a judicious cut 
in I magnitude it is possible to isolate a bright early-type galaxy sample. 

We divide the data into nine {Az = 0.1) bins (the lowest redshift z = 0.1 ± 0.05 bin contains 
very few galaxies so is discarded here). In addition to a morphological dependence, there has 
recently been quite considerable evidence in the literature for a luminosity dependence: luminous 
galaxies clustering more strongly than their fainter counterparts. Such a luminosity dependence 
has been reported by Benoist et al. (1996); Guzzo et al. (1997); Willmer et al. (1998); Small et al. 
(1999); Norberg et al. (2001); Firth et al. (2002) and Zehavi et al. (2002). The evidence for a 
dramatic increase in clustering strength with increasing luminosity for L > L* galaxies is the most 
compelling. Indeed, Norberg et al. (2002) find a factor of 2.5 difference between the clustering 
strength of L* and 4/y* galaxies and suggest that it is in fact luminosity and not type, which is the 
dominant factor (Norberg et al. measure a clustering strength for early types about 50% higher 
than that for late types at all luminosities). 

In an effort to eliminate any complicating effect of a luminosity-dependent component to uj{9), 
we decided to make a further restrictive cut to the early-type galaxies we allowed into our sample 
at each redshift. We chose to exclude all galaxies with absolute magnitude fainter than M = M^ + 1 
from our analysis. This ensures that the remaining sample has an effective luminosity Lgff ~ T* 
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(Lgff = 0.98-L* if one assumes the 2dFGRS value of L* = —19.61 (Folkes et al. 1999, See also 
Madgwick et al. 2002)). This assumes that L* for early-type galaxies does not evolve between 
z = 0 and 1. Based on our knowledge about early-type evolution with redshift this does not seem a 
grossly inaccurate assumption e.g., Lilly et al. (1995) found that the red (redder than present-day 
Sbc and hence early type) sample from the Canada-France redshift survey, was consistent with no 
change in L* between z ~ 0.8 and z ~ 0.3 (their red sample was also consistent with a change of 
at most a few tenths of a magnitude) 

Table 9 shows the number of galaxies in each redshift interval which meet our criteria. In 
estimating we once more assumed that was well described by a power-law with slope 
6 = 0.8 (=7 — 1). Slightly steeper slopes have been claimed for the spatial correlation function 
slope 7 for early-type galaxies i.e., 1.87 ± 0.07 (Loveday et al. 1995), 2.05 ± 0.1 (Guzzo et al. 1997), 
1.91 ±0.06 (Willmer et al. 1998), 1.91 ±0.06 (Shepherd et al. 2001), and 1.87±0.09 (Norberg et al. 
2002). A slightly steeper slope would have negligible effect on our conclusions, and is far from being 
the main source of uncertainty. Note also that by measuring the correlation function amplitude at 
a fixed angular scale of 1' we are measuring at increasing projected radius with redshift, ranging 
from about 140 h~^ kpc at z = 0.2 to about 330 h~^ kpc at z = 0.9. This analysis, therefore, 
assumes that any bias is not a function of scale (Magliocchetti et al. 2000). 

Figure 11 shows the logarithm of the amplitude of u){9) at 1', as a function of redshift for 
early types. The measured values of the two-point correlation function amplitude at each redshift 
may be found in Table 9. For comparison, superimposed on Figure 11 are the predictions assuming 
a correlation length of ro = 5.7 h~^ Mpc, the value determined by the 2dFGRS for local early 
types with L = A: 0.5 (Norberg et al. 2002, their Table 2). The various lines indicate the 
predictions assuming evolutionary parameter e = —1.2 (solid), 0.0 (dashed), 1.0 (dot-dashed) or 
2.0 (dotted). The resulting predictions assuming the 2dFGRS value of tq match our measurements 
of the clustering amplitude reasonably well over the whole range of redshift. Clustering fixed 
in comoving coordinates, e = —1.2, appears somewhat to overestimate the measured clustering 
amplitude. At the other extreme, rapid growth in clustering with time, e = 2.0, appears rather to 
underestimate the actual data. 

Figure 12 shows the same data as Figure 11 but plotted on a log-log scale. Here we assume a 
correlation function slope of 7 = 1.8 and a flat lambda cosmology as before. In this case, however, 
at each redshift we subtract the Gamma function and cosmological components of the two-point 
function (equations 7 and 9) from our measured value of A^. The remaining tq- and e-dependent 
portion of the correlation function amplitude [7 log ro — (3 ± e — 7 ) log(l ± z)] at each redshift is then 
plotted. The redshift baseline is too small and the uncertainties too large to solve meaningfully for 
ro and e simultaneously. We therefore solved solely for vq. The best-fit physical clustering length 
ro(€) was found to be 4.02 ± 0.22 (-1.2), 5.25 ± 0.28 (0.0), 6.55 ± 0.36 (1.0), and 8.17 ± 0.45 (2.0) 
as shown in Table 10. The lines superimposed onto Figure 12 show these best-fit estimates (solid 
for € = —1.2, dashed for e = 0, dot-dashed for e = 1, and dotted for e = 2) 
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Fig. 12 .— The logarithm of the cosmology-independent component of the angular correlation 
function uj{ 6 ) at V for L* early-type galaxies assuming 7 = 1.8. The lines show best-fit physical 
correlation length ro, for evolutionary parameters e = —1.2 (solid), 0.0 (dashed), 1.0 (dot-dashed) 
or 2.0 (dotted). Also shown are the values measured by Loveday et ah (diamond), Guzzo et ah 
(star), Willmer et ah (circle), and Norberg et al. (triangle). 
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Stronger constraints may be placed on the evolutionary parameter e if we consider our data in 
conjunction with local measurements. The Norberg et al. measurement of 5.7±0.6 (M*—0.5 <M < 
M* + 0.5) is that which most closely resembles our own study with regards to the range of absolute 
luminosity of galaxies included in the sample (M* — 0.5 < M < M* + 0.5 c.f. M < M* + 1.0 in 
this analysis). However, we also now compare with the clustering length measured locally by other 
groups for early-type galaxies (the values in parentheses show the range in absolute magnitude). 
Shepherd et al. (2001) obtained 5.45 ± 0.28 (M < M* — 1.0) from the CNOC2 survey. However, 
other groups have obtained slightly larger correlation lengths. Loveday et al. (1995) measured 
6.4 ± 0.7 (M* — 0.5 < M < M* + 0.5) for the Stromlo-APM survey, Willmer et al. (1998) measured 
6.06 ± 0.39 (M < M* -|- 5.5) for the SSRS2 survey and Guzzo et al. (1997) measured 8.35 ± 0.75 for 
the Perseus-Pisces survey. Note that the latter (Guzzo et al.) result is for a survey known to contain 
a high abundance of local clusters and therefore the high correlation length is unsurprising. Also 
plotted onto Figure 12 are the values measured by Loveday et al. (diamond), Guzzo et al. (star), 
Willmer et al. (circle) and Norberg et al. (triangle). Even with the above-mentioned caveats, it 
would appear from Figure 12 that in order to match the “local” value for early-type galaxies, the 
correlation length of such galaxies would be required to grow with time (e > 0.0). 

Equation (6) is a purely empirical model and is not necessarily expected to be valid at all 
redshifts. In particular, this model predicts monotonic evolution in tq with redshift, which is at 
odds with predictions (although these provide rather weak constraints at present) from N-body 
simulations. For example, from their semi-analytic model Baugh et al. (1999, their Figure 2) 
suggest that e = 0 is a good descriptor of the behavior of the galaxy correlation function between 
z = 0 and z = 1 i.e., clustering does not evolve in physical coordinates and there is a monotonic 
decrease in the physical correlation length with lookback time during that epoch. However, between 
z = 1 and z = 2, a minimum in the correlation length is reached. At higher redshift, Baugh et al. 
predict an upturn in the correlation length, and thus equation (6) would not describe galaxy 
clustering well in the higher redshift regime. Similar conclusions were reached from simulations 
by Kauffmann et al. (1999). Such predictions of an upturn in the correlation length at higher 
redshift are supported by the z ~ 3 measurements of Lyman break galaxies (Adelberger et al. 1998; 
Giavalisco et al. 1998) Thus, while the e model does appear to be in qualitative agreement with our 
measurements of early type clustering at z < 1, any evolution in clustering at higher redshift may 
well be poorly described by the e formalism. For now, based on our measurements, the predictions 
from semi-analytic modelling, and the requirement to match the clustering lengths measured locally 
for early-type galaxies, we assume the “null hypothesis”, namely, that clustering is fixed in physical 
coordinates. We conclude, for now at least, that the e model with e = 0.0, and a present-day 
correlation length tq = 5.25 ±0.28 provides a satisfactory fit to observations over the redshift range 
0 < z < 0.9. (See Hogg et al. (2000) for a similar conclusion of minimal evolution relative to 
stable clustering in the redshift range 0 < z < 1 although in that case comparing heterogeneous 
samples of galaxies with respect to absolute magnitude, redshift, morphological type, correlation 
function slope 7 , and cosmology). Note that for a correlation function of the form ^(r) = {r/ro)~"^, 
clustering density falls with radius as r~'^ and the cosmic mean rises with redshift as (iTz)^, so the 
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stable clustering prediction is for the proper correlation length to decrease with increasing redshift 
as ro{z) = ro( 0 )(l + . 

Interestingly, we note that the clustering strength we measure is significantly weaker than 
that found for extremely red objects (EROs) at 2 : > 1. Daddi et al. (2000, 2001) used K and 
R band imaging, McCarthy et al. (2001) used H and I band imaging, and Firth et al. (2002) 
used R and H band imaging all from the Las Campanas Infrared Survey, to select EROs (which 
are thought to be predominantly elliptical galaxies in the redshift range 1 < z < 1.5). Daddi 
et al. measured a comoving correlation length of 12 ± 3 h~^ Mpc for an effective redshift 2 : ~ 1.2. 
McCarthy et al. measured a slightly smaller comoving correlation length of 9.5 ± 0.5 h~^ Mpc. 
Using our preferred e = 0, these comoving values can be converted to physical correlation lengths 
via xo = ro(l + 2 ;)“ 0 +e- 7)/7 = ^ 0 ( 2 : = 0) which translates to ro ~ 12 — 18 h~^ Mpc. Even 

if one assumed that these galaxies underwent no evolution in comoving clustering (e = —1.2 i.e., 
xq = Xz) this would result in a present-day correlation length of ro ~ 9 —12 h~^ Mpc. Even allowing 
for these two samples containing brighter galaxies than ours it would be very difficult to reconcile 
such large correlation lengths for the ERO’s with our measurements for 0 < 2 : < 1 early types. 
Indeed, it would be difficult to reconcile such high values of correlation length with the 8.33 ± 1.82 
found locally for the brightest interval in the Norberg et al. (2002, (their Table 2)) sample. Firth 
et al. measure ro = 7.7 ±2.4 for e = —1.2 or ro = 12.1 assuming e = 0. It is possible that the width 
of the redshift distributions N{z) (derived from photometric redshifts) used in the Las Campanas 
estimates have been overestimated leading to an increase in the estimates of ro. Alternatively, the 
larger inferred values of ro may be the result of an increasing biasing with redshift (Mo & White 
1996). A yet further possibility is that ERO’s are not truly held early types but may show stronger 
clustering strengths because they are cluster ellipticals in the process of forming (Moustakas & 
Somerville 2002). Further data will be required to resolve this controversy. 


5. CONCLUSIONS 

We investigated the dependence of the two-point galaxy angular correlation function on median 
magnitude, V — I color and morphology. We found u;( 6 *) to be consistent with a power-law of slope 
—0.8 for both I and V passbands, down to our faintest limits of I = 24 and V = 25. We found 
the amplitude of co(ft} at 1', to decrease monotonically with increasingly faint median magnitude. 
We used spectroscopic redshifts measured from Cowie’s SSA22 held sample to model the galaxy 
redshift distribution as a function of apparent magnitude. We compared the measured values of 
clustering amplitude with the values predicted in an Dmo = 0.3, Qxq = 0.7 cosmology. We 
found that simple redshift-dependent models with evolutionary parameter e were inadequate to 
describe the evolution of clustering. We concluded that allowance must be made for the increasing 
proportion of later type and fainter galaxies (with weaker correlation strengths) entering our sample 
at fainter magnitudes. 

We also found a strong clustering dependence on color. Extremely blue galaxies (U — I ~ 0.5) 
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were found to have a clustering amplitude about 15 — 20 times as high as the full field population. 
This is most likely to be because many of these blue galaxies are situated at very similar relatively 
low redshift, and therefore ui{6) is minimally diluted by projection effects. Extremely red galaxies 
— / ~ 3) were found to have an clustering amplitude about 10 times as high as the full sample. 
We similarly interpreted the stronger clustering amplitude for redder galaxies to be mainly due to 
these galaxies occupying a narrow range in redshift at z = 1. The stronger signal is also due in 
part to these being early-type galaxies which locally cluster more strongly than later types. 

We then presented the first attempt to investigate redshift evolution utilizing a population of 
galaxies of the same absolute luminosity and morphological type. We used V — I color to isolate 
a sample of early-type galaxies and investigated the evolution in their clustering. By making an 
identical cut in absolute magnitude to our early-type sample at each redshift, we determined uj{9) 
for galaxies with effective luminosity Leff — L* (assuming no evolution in the luminosity function 
with redshift) in eight redshift intervals spanning z = 0.2 — 0.9. Although uncertainties were large, 
we found the evolution in the clustering of these galaxies to be consistent with stable clustering 
(e = 0). We found L* early-type galaxies to have correlation length ro = 5.25 ± 0.28 h~^ Mpc 
(assuming e = 0), a slightly higher correlation length than has been found for the local full field 
population. Our measured value of tq is in good agreement with the 2dFGRS measurement of 
correlation strength for L* early-type galaxies in the local universe. 

Over the last few years it has become increasingly apparent that galaxy clustering has a 
bivariate dependence on both morphological type and intrinsic luminosity. Clearly, if there are 
differences between the clustering of various different samples of galaxies we can immediately infer 
that at least one of the galaxy samples is a biased tracer of the underlying mass distribution. This 
paper presented the first attempt to separate the relative contributions of luminosity and type 
and to investigate the evolution in clustering of a single galaxy population. In the future, huge 
quantities of new data from galaxy redshift and large-area imaging surveys currently in progress 
will allow galaxy samples to be selected more precisely by luminosity and type. The 2dFGRS 
and SDSS are in the process of measuring spectroscopic redshifts for millions of galaxies and will 
determine, with incredible accuracy, the “local” correlation function (both slope and amplitude) as 
a function of galaxy morphological type and absolute luminosity. Preliminary findings have already 
been reported by Norberg et al. (2001, 2002); Connolly et ah (2002); Infante et al. (2002); Zehavi 
et al. (2002); Dodelson et al. (2002); Tegmark et al. (2002). At higher redshift, next generation 
redshift surveys such as DEEP2 (Davis & Faber 1998; Coil et al. 2001) or VIRMOS (Le Fevre et al. 
2001 ) will provide tens of thousands of galaxy redshifts. 

In the more immediate future, large multi-passband surveys such as the Deep Lens Survey 
(http://dls.bell-labs.com/) are measuring tens of millions of galaxies over tens of square degrees 
which will result in much more accurate photometric redshift determinations than possible in this 
study. The greater range of absolute luminosity then available (limited here to M ~ M* ± 1) will 
allow any luminosity dependence to the early-type galaxy correlation function to be determined 
more precisely as a function of redshift. Moreover, increased numbers of early-type galaxies will 
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greatly reduce uncertainties in the measurement of the amplitude and slope of the correlation 
function. Finally, the availability of more than two passbands will also allow photometric redshifts 
for late-type galaxies to be determined and a similar investigation to be undertaken into their 
clustering evolution. 
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Table 1. Field Centers and Seeing 


Field 

Pointing 

RA (J2000) 

DEC (J2000) 

1 

b 

FWHM(I) 

FWHM(V) 

Lockman 

1 

10:52:43.0 

57:28:48.0 

149.28 

53.15 

O 

bo 

CO 

0".85 


2 

10:56:43.0 

58:28:48.0 

147.47 

52.83 

0".84 

0".86 

Croth 

1 

14:16:46.0 

52:30:12.0 

96.60 

60.04 

O 

00 

o 

0".93 


3 

14:09:00.0 

51:30:00.0 

97.19 

61.57 

o 

o 

0".85 

1650 

1 

16:51:49.0 

34:55:02.0 

57.37 

38.67 

0".82 

0".85 


3 

16:56:00.0 

35:45:00.0 

58.58 

37.95 

0".85 

0".72 


Table 2. Differential Cousins I-Band Counts (N deg ^ mag 


Magnitude 

logio(N) 

'^high 

<7low 

16.0 

1.590 

0.037 

0.040 

16.5 

1.617 

0.139 

0.205 

17.0 

2.024 

0.208 

0.414 

17.5 

2.250 

0.148 

0.227 

18.0 

2.590 

0.106 

0.141 

18.5 

2.817 

0.088 

0.111 

19.0 

3.093 

0.114 

0.154 

19.5 

3.317 

0.063 

0.074 

20.0 

3.541 

0.044 

0.049 

20.5 

3.737 

0.020 

0.021 

21.0 

3.905 

0.025 

0.027 

21.5 

4.078 

0.020 

0.021 

22.0 

4.236 

0.032 

0.035 

22.5 

4.395 

0.037 

0.041 

23.0 

4.551 

0.026 

0.027 

23.5 

4.696 

0.020 

0.021 

24.0 

4.797 

0.024 

0.025 











Table 3. Best Fits to Counts 


Passband 

Range 

Slope 

Intersection 

I 

[21.0-23.5] 

0.31 ±0.01 

-2.68 ±0.25 

V 

[22.0-24.5] 

0.43 ± 0.02 

-5.83 ±0.45 


Table 4. Comparison of Number Count Slopes 


Passband 

Range 

Slope 

Source 

I 

[21.0-23.5] 

0.31 ±0.01 

This Work 

I 

[20.0-24.0] 

0.35 ±0.02 

McCracken et al. (2001) 

I 

[21.0-25.0] 

0.33 

Metcalfe et al. (2001) 

I 

[22.5-25.5] 

0.31 ±0.02 

Arnouts et al. (1999b) 

I 

[19.0-22.5] 

0.34 ± 0.03 

Driver et al. (1994) 

V 

[22.0-24.5] 

0.43 ± 0.02 

This Work 

V 

[22.0-24.5] 

0.43 ± 0.03 

Woods & Fahlman (1997) [averaged] 

V 

[22.0-24.5] 

0.404 ± 0.015 

Smail et al. (1995) 

V 

[20.5-23.0] 

0.41 ±0.01 

Driver et al. (1994) 
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Table 5. Differential V-Band Counts (N deg ^ mag 


Magnitude 

logio(N) 

'^high 

O'low 

17.5 

1.667 

0.178 

0.306 

18.0 

1.807 

0.138 

0.204 

18.5 

2.037 

0.240 

0.580 

19.0 

2.309 

0.083 

0.103 

19.5 

2.606 

0.108 

0.144 

20.0 

2.843 

0.097 

0.125 

20.5 

3.035 

0.095 

0.122 

21.0 

3.255 

0.065 

0.077 

21.5 

3.422 

0.057 

0.066 

22.0 

3.644 

0.049 

0.055 

22.5 

3.849 

0.037 

0.040 

23.0 

4.059 

0.035 

0.038 

23.5 

4.298 

0.027 

0.029 

24.0 

4.521 

0.031 

0.034 

24.5 

4.701 

0.030 

0.032 

25.0 

4.796 

0.040 

0.044 


Table 6. The Dependence of ^^(1^) on Median Magnitude 


Range 

Median 

"^gal 

logio ^a;(l') 

20.0 <I < 21.0 

20.60 

5235 

-1.21 ±0.08 

21.0 <I < 22.0 

21.60 

11535 

-1.49 ±0.04 

22.0 <I < 23.0 

22.59 

23842 

-1.80 ±0.04 

23.0 <I < 24.0 

23.57 

47141 

-2.19 ±0.07 

21.0 < D < 22.0 

21.62 

2598 

-1.13 ±0.06 

22.0 < D < 23.0 

22.61 

6883 

-1.30 ±0.05 

23.0 < D < 24.0 

23.63 

19411 

-1.72 ±0.05 

24.0 < D < 25.0 

24.58 

47057 

-2.26 ±0.09 
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Table 7. Best-fit Redshift Scale Parameter zq as a Function of Magnitude 


Range 

Raw 

Corrected 

20.0 <I < 20.5 

0.15 

0.15 

20.5 <I < 21.0 

0.20 

0.20 

21.0 <I < 21.5 

0.21 

0.22 

21.5 <I < 22.0 

0.22 

0.26 

22.0 <I < 22.5 

0.26 

0.36 

22.5 <I < 23.0 

0.30 

0.41 

23.0 <I < 23.5 

0.34 

0.43 


21.5 < P < 22.0 

0.13 

0.13 

22.0 < P < 22.5 

0.15 

0.15 

22.5 < P < 23.0 

0.17 

0.17 

23.0 < P < 23.5 

0.20 

0.20 

23.4 < P < 24.0 

0.24 

0.26 

24.0 < P < 24.5 

0.30 

0.37 


Table 8. The Dependence of onV — I Color for Galaxies in the Range 20.0 < I < 23.0 


Color Range 

-^gal 

logio 

0.0 < P 

- / < 0.5 

136 

0.06 ± 0.45 

0.5 < P 

- / < 1.0 

4209 

-1.27 ±0.08 

1.0 < P 

- / < 1.5 

13511 

-1.64 ±0.08 

1.5 < P 

-I < 2.0 

13587 

-1.50 ±0.04 

2.0 < P 

-I <2.5 

5757 

-1.20 ±0.08 

2.5 < P 

- / < 3.0 

4211 

-1.03 ±0.12 

3.0 < P 

-I <3.5 

1833 

-0.78 ±0.08 
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Table 9. A^{1') for Early-Type Galaxies as a Function of Redshift 


Redshift 

^gal 

logio ^a;(lO 

0.2 ±0.05 

136 

0.05 ±0.36 

0.3 ±0.05 

366 

-0.21 ±0.07 

0.4 ±0.05 

569 

-0.54 ±0.11 

0.5 ±0.05 

559 

-0.53 ±0.15 

0.6 ±0.05 

389 

-0.60 ± 0.26 

0.7 ±0.05 

551 

-0.38 ±0.10 

0.8 ±0.05 

575 

-0.47 ±0.16 

0.9 ±0.05 

237 

-0.28 ±0.12 


Table 10. 


Best-Fit Physical Correlation Length rg and Reduced for Various Values of 
Evolutionary Parameter e. 


e 

ro 

X^/dof 

-1.2 

4.02 ±0.22 

1.01 

0.0 

5.25 ±0.29 

1.59 

1.0 

6.55 ±0.36 

2.70 

2.0 

8.17 ±0.45 

4.40 



